Competing orders in one dimensional spin 3/2 fermionic systems 
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Novel competing orders are found in spin 3/2 cold atomic systems in one-dimensional optical traps 
and lattices. In particular, the quartetting phase, a four-fermion counterpart of Cooper pairing, 
exists in a large portion of the phase diagram. The transition between the quartetting and singlet 
Cooper pairing phases is controlled by an Ising symmetry breaking effect in one of the spin channels. 
The singlet Cooper pairing phase also survives in the purely repulsive interaction regime. In addition, 
various charge and bond ordered phases are identified at commensurate fillings in lattice systems. 
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Optical traps and lattices open up a whole new direc- 
tion in the study of strongly correlated systems in cold 
atomic physics. In particular, they provide a controllable 
way to investigate high spin physics by using atoms with 
hyperfine spin multiplcts. For example, polar and ferro- 
magnetic condensations in spin-1 bosonic systems (e.g. 
23 Na and 87 Rb) have been extensively studied both ex- 
perimentally 0, Q and theoretically Q. Spin nematic 
orders and the exotic spin liquid behavior are arousing 
much interest |j, la, |g, |j| m Mott-insulating phases. On 
the other hand, high spin fermionic systems also exhibit 
many properties different from those in the usual spin 
1/2 systems. For example, Cooper pairing shows new 
structures with total spin 5 > 2 in the s-wave channel 



High spin systems also provide a wonderful opportu- 
nity to investigate the multi-particle clustering (MPC) 
instability and its competition with Cooper pairing in 
cold atomic physics. Taking into account the recent ex- 
citing achievement of fermionic superfluids by using Fes- 
hbach resonances ^(J> the MPC instability will be one of 
the next focus directions. For example, the three-body 
Efimov bound states in bosonic systems have attracted 
wide attention In fermionic systems, the MPC insta- 
bility is forbidden in two-component spin 1 /2 systems due 
to Pauli's exclusion principle, but is allowed in high spin 
systems with multiple components. Knowledge learned 
from the MPC phase in cold atomic systems will also shed 
light on the a-particle formation 12] in nuclear physics, 
which is a four-fermion bound state. Some previous theo- 
retical works have studied the MPC phase in the SU (N) 
symmetric fermionic models by the Bethe ansatz at ID 
pL3| and a variational method at high dimensions [T^| . 
However, the stability of the MPC phase remains to be 
an open problem when systems do not possess the SU (N) 
symmetry, and the nature of the transition between the 
MPC and Cooper pairing phases has not been clarified 
before. 

Spin 3/2 systems provide an ideal starting point to 
study the simplest MPC phase in fermionic systems, i.e., 
the quartetting phase, which can be realized by using 
132 Cs, 9 Be, 135 Ba, 137 Ba atoms. In this article, we in- 



vestigate the quartetting phase in one dimensional (ID) 
spin 3/2 systems where its strong coupling nature can 
be handled by applying the methods of bosonization and 
renormalization group (RG). The investigation is greatly 
facilitated by the generic 190(5) symmetry recently iden- 
tified in spin 3/2 systems |l5j, [lfj. We find rich phase 
structures, including a gapless Luttinger liquid phase and 
two distinct spin gap phases at incommensurate fillings. 
One of these spin gap phases is characterized by the quar- 
tet formation, while the other is dominated by Cooper 
pairing. The transition between them is controlled by 
an Ising symmetry breaking in one of the spin channels. 
Both quartets and Cooper pairs can further undergo ei- 
ther the quasi-long range ordered (QLRO) superfluidity 
or charge density wave (CDW) instabilities. In contrast, 
the 5/7(4) symmetry in the previous study 0] requires 
four particles to form an 5C/(4) singlet, thus in ID only 
the quartetting phase is allowed and Cooper paring is 
excluded. Furthermore, various charge and valence bond 
orders are identified at commensurate fillings. 

We begin with phase structures at incommensurate 
fillings. After linearizing the spectra around the Fermi 
wavevector kf, we decompose fermion operators into 
right and left moving parts as i/; a = ip R , a e lk,x + 
^L, a e~ lkfX {a = ±|,±i). The right (left) moving cur- 
rents are classified into 50(5) 's scalar, vector and tensor 
currents as J R{L) (z) = ip RWia (z)ip R (L), a (z), J R(L) (z) = 

I^Ll^lryHW,^)! 1 < a < 5), and J R \ L) (z) = 
¥ R (L), a ( z ) Ta ^R(L)A z )( l <a<b<5), where T a ,T ab 
are the 4x4 Dirac matrices defined in Ref. ^5|. Classi- 
fied in terms of the usual spin SU (2) group, the scalars 
Jr(l) are charge currents, the 5-vectors J R r L \ are spin- 
nematic currents with spin 5 = 2, and the 10-tensors 
Jr(l) con tain two degenerate parts with spin 5=1,3. 

Spin 3/2 systems are characterized by two independent 
s-wave scattering parameters go and 172 in the total spin 
singlet (St — 0) and quintet (St = 2) channels, respec- 
tively. Taking these into account, the low energy effective 
Hamiltonian density reads 

Ho = v f {^J R jR + ^(J%J% + jgjg) + (R^L)}, 
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'Hint — ^tJrJl + 9vJrJl + 9tJR 'JL b 



(1) 



where g c ,v t are effective dimensionless coupling constants 
with their bare values given by 

2g c = 9o + 5.g 2 , 2g v = g Q - 3g 2 , 2g t = -(g + g 2 ). (2) 

We neglect the chiral interaction terms. Eq. \T\is closely 
related to the extensively studied two-coupled spin 1/2 
Luttinger liquids [1J], UJI Uj| • In the latter case, Fermi 
wavevectors in two bands are usually different, while they 
are the same for all of the four components in the former 
case, thus much more low energy inter-band interactions 
are allowed. An even larger symmetry SU (4) can be 
obtained by fine tuning go = g 2 , i.e., g v = g t , which 
means that the four spin components are equivalent. In 
this case, and J R \ L \ form the generators of the 

SU(4) group. Then Eq. U reduces to the single chain 
SU(A) spin-orbit model [ffi . 

The RG equations can be derived through the current 
algebra and operator product expansion techniques [2~i| . 
The charge current satisfies the U(l) Kac- Moody algebra, 
and thus g c is not renormalized at one loop level due to 
the absence of Umklapp terms. The vector and tensor 
currents in spin channels form the SU (4) Kac-Moody 
algebra, which gives rise to 



dg v _ _4_ dg t _ J_ 

d\n(L/a) ~ 2ir 9v9u d\n(L/a) ~ 2ir 



(3 5 2 + «? 2 ), (3) 



where L is the length scale and a is the short distance 
cutoff. The SU (4) symmetry is preserved in the RG pro- 
cess along the line g v = gt- Eq. [3] can be integrated as 
\9t ~ 9v\ = c \9v | 3 ^ 2 (c ; constant) with the RG flows as 
shown in Fig^ The parameter space (g v ,gt) is divided 
into three phases as A) the gapless Luttinger liquid, B) 
and C) two different spin gap phases with the formations 
of quartets and singlet pairs respectively, which will be 
clarified below. Phase A lies between g t — g v < (line 1) 
and gt = —g v < (line 2), where RG flows evolve to the 
fixed point (FP) g v = gt = 0. Phase B is bounded by line 
2 and g v = (line 3) where RG flows evolve to the SU (4) 
marginally relevant FP of g v — g s — > +00 (line 4) . Phase 
C is symmetric to B under the reflection g v — > —g v . RG 
flows evolve to the line 5 with —g v — g t — > +oo. The 
boundary between B and C, i.e., line 3, is controlled by 
the unstable FP of g v = 0, gt — > +oo. From the relations 
in Eq. [21 we replot the phase diagram in terms of the s- 
wave scattering lengths go and g 2 as shown in Fig. [21 The 
boundary between phase A and C, i.e., go = 92 > 0, is ex- 
act from the SU(A) symmetry regardless of the one-loop 
approximation. On the contrary, boundaries between A 
and B, B and C will be modified in higher order pertur- 
bations. 

To clarify the nature of each phase, we employ 
the Abelian bosonization. The bosonization iden- 
tity reads ipR( L ) a (x) = g a /V2ira exp{±i v /7r(0 C[ (a;) ± 
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FIG. 1: RG flows in the spin channels. Various phases are 
determined as: A) the gapless Luttinger liquid phase; B) the 
quartetting phase with QLRO superfluidity at K c > 2 or 2k /- 
CDW at K c < 2; C) the singlet pairing phase with QLRO 
superfluidity at K c < 1/2 or 4fc/-CDW at K c > 1/2. They 
are controlled by the FPS of (0, 0),(+oo, +oo) (line 4), and 
(— oo,+oo) (line 5) respectively. Phase boundaries (line 1, 2, 
3) are marked with dashed lines. 



9 a (x))}(a — ±|,±i) where 7]s are the Klein factors. 
Boson fields <f> a and their dual fields 9 a are conve- 
niently reorganized into 4> C (9 C ) in the charge channel, 
and 4'v(9 v ),(f>ti{9ti),4 > t2{0t2) in the spin channels via 

<pc,v = (0| ± <t>\ ± <\>-\ + 0-§)/2,0tl, f 2 = (</>§ T 4>\ ± 

<p_i — cj)_3 )/2. Similar expressions also hold for 9s. The 
quadratic part of the Hamiltonian density is standard 



(v = c,v,h,t 2 ) 



with Luttinger parameters K v and velocities v v in each 
channel. The non-quadratic terms are summarized as 

Hint = -777— T^l cos V47T0tl + COS V^TT(f> t2 } 

2{na) z 



x { (9t + 9v) cos v4ar(j> v + (g t - g v ) cos v47T^t, } 

gt 



2(iray 



COS V 4:lT(j>tl COS V 47T0 t2 



(•5) 



with the convention of Klein factors as 77377177-177-3 =1. 

Various order parameters are needed to characterize 
phase structures, including the 2fc/-CDW operator N, 
the spin density wave (SDW) operators in the SO(5) vec- 
tor channel A^ a and their tensor channel version N ab , the 
singlet pairing operator 77 and its quintet counterpart x°, 
the 4&/-CDW operator 04k f , c dw, as well as the quartet- 
ting operator O q t- They are defined as 



3 
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FIG. 2: Phase diagram in terms of go and <?2- We also show 
the charge and bond ordered states at quarter filling in lattice 
systems. Phase A is the SU(4) gapless spin liquid. Both phase 
B and C split into two parts. With g u — > +00, B.l) quartets 
with both bond and charge orders, C.l) dimerization of spin 
Peierls order. With g u —* —00, B.2) CDW phases of quartets 
and C.2) CDW of singlet Cooper pairs. Boundaries among 

A, B and C are marked with solid lines, and those between 

B. l and B.2, and C.l and C.2 are sketched with dashed lines. 



4k f ,cdw 



quar 



(6) 



where R = T\T^ |l9j is the charge conjugation matrix 
and e a f3js is the rank-4 Levi-Civita symbol. 

Now we are ready to discuss the nature in each phase. 
In the Luttinger liquid phase A, all the cosine terms are 
marginally irrelevant. As shown in Fig. [21 the Luttinger 
phase exists in the repulsive region with go > g 2 > 0, thus 
K c < 1. As a result, only 2fc/-CDW and SDW suscep- 
tibilities diverge with scaling dimensions Ajy = Ajy^ = 
A N ab = (if c + 3)/4. Taking into account the chiral inter- 
action terms, two different spin velocities exist in the vec- 
tor and tensor channels as exhibited in the one-particle 
correlation functions (\i/>R( L ), a {%, r )^fi(L) ia! ( ' °)D 



1 



1 / a y° , 7 - 

Fix) 3 \ \v r T - ix\J 



(v c T^fix)i (v v t =f ix)* (v t T =f ix)i M^c 

with 7c = \(K c + l/K c -2). 

Phase B is characterized by spin gaps in <f> v , <f>ti and 
4>t2 channels. Its SU (4) FP of g v = gt — > +00 means the 
formation of the SU(4) invariant quartets as follows. In 
the ground state, cosine operators acquire non- vanishing 
expectation values as (cos V^br^u) = (cos %/47n/> t i) = 
(cos V^TT(f>t2), whose classical values are just 1 or —1 as 
related by a Z 2 symmetry 22]. We fix the gauge by 
choosing (<f> v ) = (<j>tx) = {4>t2) — 0. By checking scal- 
ing dimensions of various order parameters in Eq. EJ we 
conclude that the competing instabilities are the QLRO 
superfluidity O quar and CDW of quartets. Because the 
average distance between two quartets is d = Tr/kf, this 



CDW is of 2k f type. Their expressions are reduced to 
O qU ar oc e 2l ^ e " and N cx e 1 ^^" respectively, with scal- 
ing dimensions A quar — 1/K C and An — \K C . The 
leading instability is O quar at K c > 2, and N at K c < 2 
respectively. We extend the quartets formation regime 
from the previous Bethe-ansatz results along SU (4) line 
4 0| to the entire phase B. On the other hand, correla- 
tions of pairing operators 77^ and decay exponentially. 

Phase C is controlled by the FP of — g v — g t — ► +00 
which is characterized by the singlet Cooper pairing, and 
possesses another SU{4) symmetry denoted as St/' (4). 
Its right (left) generators J' R i L \ belong to the SU (4) fun- 
damental (anti-fundamental) representations defined as 



jf,ab jab j',a ja j/,ab jab jt,a ja 

J R ~ J R > J R ~ J R> J L — J L i J L ~ ,J L- 



(8) 



It is different from the St/ (4) symmetry in phase B where 
both right and left generators are in the fundamental 
representation. The singlet Cooper pair is invariant un- 
der this St/' (4) symmetry but not the St/(4) symmetry 
in phase B. We choose the pinned bosonic fields to be 
(0ti) = (<fe) = (6v) — 0, where the dual field 9 V instead 
of 4> v is pinned as in Phase B. Again, the superfluidity 
and CDW orders of pairs compete. Because the aver- 
age distance between adjacent pairs is d = n/(2kf), this 
CDW is of the 4/c/ type. These two orders are reduced 

with scaling 
, = K c . The 



to 77+ cx e l ^ 9c and O ikj , c dw oc 
dimensions as A„ = 1/ (4AT C ) and A 



/4lT0 c 



leading instability is if at K c > 



4kf,cd 

1/2, and 4kfiCdw at 
K c < 1/2 respectively. Remarkably, due to the presence 
of the spin gap, the singlet pairing instability dominates 
even in the purely repulsive region at go > g 2 > when 
K c > 1/2 as shown in Fig. [21 This is similar to the 
situation in high T c superconductivity. 

The transition between quartetting and pairing phases 
is controlled by an Ising duality in the 9 v (<p v ) channel. 
Near the phase boundary, 0(1,2 are pined, thus the resid- 
ual interaction becomes 



H r 



A 



2(7ra) 2 
X cos 



{(gt + 9v) cos VAtt<P v + (g t - g v ) 



(9) 



where A = —^((cosy/Anftti) + (cos V^n(f> t2 )). The 
singlet pairing operator cx ip^tps — ip\ip-i cx 

2 ~2~ 2 ~2~ 

e -iV^o e ^ e -i^e v + e iV*o v ^ thug s [^ x q v is t he relative 

phase between the pairs of A{ = t/j'sV'I.s and A\ = 

22 

V'iV'Ii > an( i is the vortex field dual to 9 V . The tran- 

2 2 

sition can be viewed as the phase locking problem in a 
two-component superfluids. In the pairing phase, 9 V can 
be pinned at either or ^JH by the internal Josephson 
term cos \/~Ak9 v . This symmetry is Z 2 , and thus the pair- 
ing phase with pinned 9 V is in the Ising ordered phase. 
In contrast, the quartetting phase is the Ising disordered 



4 



phase where the dual field 4> v is locked instead. The Ising 
nature of the transition is also clear by representing Eq. 
[5] in terms of two Majorana fermions [23j as 



a=l 



£ is always off-critical as gt —* +00, while g v is the mass 
of £ 2 which changes sign across the boundary. Thus the 
quartetting and pairing phases are Ising disordered and 
ordered phases for £ 2 , respectively. 

Next we consider the charge and valence bond ordered 
state as depicted in Fig. |21 at quarter-filling, i.e, one 
particle per site. The 8k t Umklapp term appears, and 
decouples with spin channels, which can be bosonized as 



TL 1 



2(7ra) 2 



8k s x), (10) 



where g u is at the order of 0{g$, #f) a * the bare level. 
At K c < 1/2, this term is relevant and opens the charge 
gap. The real parts of N and Oik t ,cdw describe the usual 
CDW orders, while their imaginary parts mean the bond 
orders, i.e., the 2k f and 4k t spin Peierls orders. Phase 
A remains gapless in spin channels as described by the 
SU(4) Heisenberg model. The quartetting phase B splits 
into two parts B.l and B.2 with g u — > +oo(— 00) re- 
spectively. In B.l, Im04k f ,cdw, ReN and ImN are long 
range ordered. The ££7(4) singlet quartets exhibit both 
charge and spin Peierls orders. Instead, the 2fc/-CDW of 
quartets becomes long range ordered in B.2. Similarly, 
the singlet pairing phase C splits into C.l and C.2 as 
g u —>■ +oo(— 00) respectively. In C.l, Im04k f , c dw, i.e., 
the spin Peierls order is stabilized. Instead, the CDW 
of singlet pairs becomes long range ordered in C.2. The 
boundaries between phase B.l and B.2, phase C.l and 
C.2 are determined by the bare value of g u o = as 
sketched in Fig. However, due to the non-universal 
relation between g U Q and go. 2 , the exact boundaries are 
hard to determine. 

The model discussed above can be realized by load- 
ing spin 3/2 fermions into one-dimensional optical tubes. 
The quartetting phase can be probed by the radio- 
frequency spectroscopy to measure the excitation gaps of 
the successive quartet breaking process as 4 — > 1 + 3 — > 
1 + 1 + 2— > 1 + 1 + 1 + 1. In contrast, in the pairing 
phase, only one pairing breaking pro cess of 2 — > 1 + 1 
exits as measured experimentally [24J. The charge and 
bond ordered states break translational symmetry. The 
periodicity of phases CI and C2 is two lattice constants, 
while that of Bl and B2 is doubled. These can be de- 
tected by using the elastic Bragg spectroscopy [j25| . 

In summary, we have investigated the global phase dia- 
gram in ID spin 3/2 systems, including the gapless Lut- 
tinger liquid, the spin gapped quartetting and Copper 
pairing phases. The competition between the quartet- 
ting and pairing phases is found to be controlled by an 



Ising duality. Both quartets and pairs can either undergo 
QLRO superfluidity or CDW instabilities, depending on 
the value of Luttinger parameter K c . The QLRO sin- 
glet pair superfluidity can survive in the purely repul- 
sive interaction regime. The Mott-insulating phases at 
commensurate fillings exhibit various charge and valence 
bond orders. 



Note added After the paper was submitted for pub- 
lication, we became aware of the work by Lecheminant, 
Azaria, and Boulat [2f|, where the two spin gap phases, 
and Mott phases at quarter filling are obtained indepen- 
dently. 
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